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Midge distribution patterns at multiple spatial scales 
in shallow Finnish lakes revealed through sedimentary records 


Tomi P. Luoto, Liisa Nevalainen & Veli-Pekka Salonen 


1. Introduction 


Luoto, T. P., Nevalainen, L. & Salonen, V.-P. 2012: Midge distribution patterns at 
multiple spatial scales in shallow Finnish lakes revealed through sedimentary re- 
cords. — Entomol. Fennica 23: 83—94. 


This study was based on sedimentary midge (Diptera: Nematocera) assemblages 
from multilake datasets along environmental transects from Finland (regional), 
southern Finland (semiregional), and Helsinki district (local) and an intralake 
dataset from eastern Finland (site-specific). The aim was to examine scale-de- 
pendencies in midge distribution. The results imply that distribution and abun- 
dance of midge taxa are related to scale: on the regional scale the forcing factors 
are related to prevailing climate conditions, on semiregional scale they are related 
to water quality, on more local scales predation pressure is the key variable and on 
site-specific scales habitat characteristics determine the species assemblages. Al- 
though the number of study sites was not equal and not all environmental parame- 
ters were possible to measure from all spatial scales, it is apparent that caution is 
required in midge-based environmental assessments, because changes in faunal 
composition are driven by factors operating at different spatial scales. 
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recent years, can sometimes be beneficial com- 
pared to traditional ecological monitoring. In the 


For the purposes of modern aquatic ecology, 
“palaeoecological” multilake datasets can be 
very informative, because they are able to pro- 
vide large amounts of data on the modern distri- 
bution and abundance of taxa. Fossil assemblages 
in surface sediment samples, which provide inte- 
grated species data from a time period of several 


case of invertebrates that vary seasonally in their 
occurrence, such as larval midges (Insecta: Dip- 
tera: Nematocera) and water fleas (Crustacea: 
Cladocera) (Einarsson ef al. 2002, Tatosova & 
Stuchlik 2006, Nevalainen 2008, Nevalainen & 
Sarmaja-Korjonen 2008), the use of ecological 
monitoring to establish modern distribution data 
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in lakes would be very time- and resource-con- 
suming and would also provide information on 
fewer taxa than the approach using sedimentary 
remains (Saros 2009, Nevalainen 2010). How- 
ever, the fossil taxonomy of midges suffers from 
lower resolution compared to living larvae or pu- 
pal exuviae (Luoto & Raunio 2011) and in addi- 
tion, the traditional multilake (single-samples) 
approach does not distinguish which species in- 
habit particular habitats within a lake basin. Re- 
cent results from multisample studies within a 
single lake have shown that it is possible to distin- 
guish habitat-specific information from surface 
sediment fossil data, as it appears that the sedi- 
mentary midge assemblages represent fauna that 
lived in close proximity to where their remains 
were deposited after moulting (instars) or death 
(Luoto 2010). There has also been recent major 
advancement in the identification of fossil midge 
remains and many of the taxa can now be identi- 
fied to species-type level (Heiri et al. 2004, 
Brooks et al. 2007). 
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Fig. 1. The predominat- 
ing forcing factors on 
midge distribution in 
Finland (60-70 °N) at 
multiple spatial scales. 
The lakes (regional 
dataset) are marked in 
open circles and sites 
(site-specific dataset) in 
black dots. The lakes in 
the semiregional and 
local datasets are all in- 
cluded in the regional 
dataset. 


Macroinvertebrate communities respond to 
multiple environmental factors, many of which 
are scale-related (Vinson & Hawkins 1998, Puntí 
etal. 2009). However, the ecological responses to 
scale are usually not linear (Krummel et al. 1987) 
and interpretation across scales have several 
problems, such as different mechanisms operat- 
ing at different spatial scales (McIntyre 1997) and 
regional and local effects may be differing (Roth 
etal. 1996). Due to these problems, consideration 
of scale can affect estimation of community struc- 
ture (Fuisz & Moskat 1992) and, therefore, must 
not be ignored in ecology or palaeoecology. Al- 
though it is well known that multiple factors oper- 
ating at various spatial and temporal scales shape 
species distribution (Levin 1992), less is known 
about how the determinants of the distributions 
(environmental forcing factors) vary across spa- 
tial scales, especially in the case of insects (Hortal 
et al. 2010). Importantly, changes in faunal com- 
position driven by factors operating at different 
spatial scales can cause serious problems in 
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Table 1. Mean values and range (in brackets) of the measured environmental variables in the regional, 
semiregional, local, and site-specific datasets and the lengths of DCA axis 1 and 2 based on the midge assem- 
blages in each of the datasets. The single-point oxygen and pH measurements in the regional dataset were 
measured from only 51 lakes. Water colour and presence of fish (fish index, Fl) were measured only from the lo- 
cal dataset and macrophytes (sample distance to littoral vegetation, DLV) and stream flow (lotic index, LI) only 


from the site-specific dataset. 








Regional Semiregional Local Site-specific 
N 82 (lakes) 30 (lakes) 9 (lakes) 34 (sites) 
Sampling depth (m) 3.3 (0.5-9.0) 3.2 (1.0-9.0) 3.7 (1.5-5.9) 3.1 (0.5-7.3) 
Mean air T, (°C) 14.9 (11.3-17.1) 16.5 (15.9-17.1) 16.6 (16.4-17.1) 15.3 
Elevation (m) 82.9 (11.3-404.0) 84.1 (11.3-134.1) 66.0 (14.6—97.3) 188.4 
Surface area (km’) 0.3 (0.003-6.94) 0.24 (0.01-6.00) 0.04 (0.01-0.14) 0.3 
Oxygen (mg I’) 5.6 (0.5—11.8) 5.0 (0.5-11.8) 8.9 (5.6—11.8) 10.9 
pH 6.2 (3.8-9.3) 6.3 (4.4-9.3) 5.5 (4.6-7.3) 6.8 
Conductivity (uS cm) 38.3 (8-312) 77 (8-312) 37 (16-128) 12 
Water colour (Pt mg I’) No data No data 55 (25-125) No data 
Presence of fish (FI) No data No data 0/1 No data 
Macrophytes (DLV m) No data No data No data 185 (0-270) 
Stream flow (LI) No data No data No data 0/1 
Length of DCA axis 1 2.9 2.0 1.7 2.3 
Length of DCA axis 2 32 1.9 1.2 3.2 





midge-based environmental assessments, for ex- 
ample in lake management at the local/site-spe- 
cific scale (reference conditions) or in climate 
change assessments at the global/regional scale 
(temperature dynamics). In addition, understand- 
ing scale-dependencies can be useful in conserva- 
tion planning of insects (Cabeza et al. 2010). 

In this study, we examine distribution patterns 
of fossil midges in shallow lakes in Finland. 
Midge-based datasets from the same region, but 
at different spatial scales, i.e. regional (Luoto 
2009a), semiregional (Luoto & Salonen 2010), 
local (Nevalainen & Luoto 2010), and site-spe- 
cific (Luoto 2010) datasets, are compared for 
their most important environmental forcing fac- 
tors on midge taxa distribution and abundance. 
The datasets represent species-environment rela- 
tionships at two grain sizes: the lake (regional, 
semiregional, and local datasets) and site (site- 
specific dataset) levels. The hypothesis here is 
that different environmental factors have a domi- 
nant role at particular spatial scales. 


2. Study sites 


The study area is located in Finland, ranging from 
60°13’ to 69°53’ N and 22°00’ to 30°13’ E (Fig. 


1). The study sites span a 1,080 km latitudinal 
transect, from boreal coniferous forests in the 
south to tundra vegetation in the north. The exam- 
ined 83 lakes are small and shallow (Table 1). The 
altitude of the lakes varies from 11 to 404 m above 
sea level (a.s.1.) and the mean July air temperature 
(7,,,) ranges from 11.3 to 17.1 °C (Table 1). The 
mean annual precipitation varies within the area 
from 400 to 660 mm (Hyvärinen 1986). The sur- 
face sediment dataset representing the regional 
scale was collected from 82 of these sites. The 
dataset representing the semiregional scale was 
collected from 30 lakes in southern Finland 
(60°13’-63°05’ N, 22°00’-30°13’ E) and the 
scaling is based on geographical distribution of 
the lakes. These lakes constitute from the south- 
ern lakes in the regional dataset. The vegetation in 
the catchments of the 30 lakes in this dataset con- 
sists of boreal coniferous forests. Some of the 
lakes are located in urban areas. In general, the 
study area at the semiregional scale is more influ- 
enced by human activity than that examined at the 
regional scale, although many of the sites have re- 
mained close to their pristine state. The limno- 
logical parameters: pH, conductivity and dis- 
solved oxygen (DO) in the regional and 
semiregional datasets were highly variable. More 
detailed site descriptions and environmental data 
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are given in the original publications (Luoto 
2009a, Luoto & Salonen 2010). 

The nine study lakes representing the local 
scale are distributed in close proximity to each 
other in southern Finland, near the Helsinki dis- 
trict (Fig. 1). The catchments consist of boreal co- 
niferous forests and the lakes vary from oligo- to 
eutrophy. They are all basins without major inlets 
or outlets. Six of the lakes are currently inhabited 
by fish (Perca fluviatilis) and three lakes are cur- 
rently without fish populations. These lakes also 
occur in the regional and semiregional datasets. 
More detailed descriptions of the sites and envi- 
ronmental data in the local dataset are given in a 
previous publication (Nevalainen & Luoto 
2010). 

Samples of the site-specific dataset were de- 
rived from a dystrophic lake, Pieni-Kauro, situ- 
ated in Kuhmo, eastern Finland (Fig. 1). The area 
is sparsely inhabited and the catchment is charac- 
terized by a drumlin landscape and wilderness- 
like coniferous forests. The northern shore of the 
lake is covered by wetlands and the Saavanjoki 
River enters the lake from the north, while the 
main outlet is on the western side. Lake Pieni- 
Kauro is part of a larger lake complex connected 
by small rivers that originates in Russia and even- 
tually drains into the Gulf of Bothnia in the Baltic 
Sea. A detailed site description, including aquatic 
vegetation patterns, is available in a previous 
publication (Luoto 2010). The samples were se- 
lected to represent different habitat types (sub- 
strate, vegetation, stream flow) and a water depth 
gradient (maximum depth 7.3 m). In all, 34 samp- 
les were collected, of which three were riverine 
and 31 lacustrine. Lake Pieni-Kauro is not in- 
cluded in the other datasets. 


3. Material and methods 


The surface sediment samples of the regional 
(N=82), semiregional (V=30), and local (V=9) 
datasets were obtained through the lake ice with a 
Limnos gravity corer (Kansanen et al. 1991) be- 
tween February and April 2005. The surface sedi- 
ment sampling of Lake Pieni-Kauro (N=31) and 
its inlet, the Saavanjoki River (N=3), (site-spe- 
cific dataset) was performed in October 2008 
from a boat, also using the Limnos corer. The sur- 
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face sediment samples (topmost 0—1 cm) used in 
this study represent the recent years of sedimenta- 
tion. The fossil faunal assemblages in these samp- 
les can be considered as analogues of modern rel- 
ative abundances (Frey 1960, Bunbury & Gajew- 
ski 2008, Nykänen et al. 2009). 

Limnological measurements were performed 
as single measurements taken in situ during the 
collection of the sediment samples with an Orion 
Model 1230 pH/mV/ORP/conductivity/dissol- 
ved oxygen/salinity/temperature meter (Thermo 
Fisher Scientific, Waltham, MA, USA) before 
sediment sampling to avoid the resulting distur- 
bance of the water column. The oxygen concen- 
tration was recorded ca. 0.5 m from the bottom at 
the sampling site and conductivity and pH from 
the epilimnion. The colour of water represents the 
platinum concentration in a platinum cobalt chlo- 
ride solution that has the same colour as the samp- 
le. The depth was measured at the sampling point 
and is in fact the water depth at the sampling site, 
which does not mean the same as the maximum 
lake depth. The multilake dataset samples were 
taken from the sublittoral areas rather than from 
the deepest point to better include the littoral taxa. 
The modern mean air 7,,, was estimated for each 
lake using a geographical information system 
(GIS)-based method. The temperature data were 
provided by the Finnish Meteorological Institute 
and were based on the 1971—2000 climate nor- 
mals, which take into account all Finnish meteo- 
rological data. The temperature data were inter- 
polated from the data of the meteorological sta- 
tions toa 10 x 10 km grid, taking into account fac- 
tors such as elevation and lake coverage. The fish 
status of the lakes in the local dataset was deter- 
mined by test fishing, diving, and available litera- 
ture (Nyberg et al. 2010). The distance from the 
littoral vegetation (DLV) was measured and de- 
termined by examining aerial photographs and 
field observations and stream flow was deter- 
mined simply by dividing the samples following 
their lotic or lentic distribution. For more details 
and other measured environmental parameters, 
consult the previous publications (Luoto 2009a, 
2010, Luoto & Helama 2010, Luoto & Salonen 
2010, Nevalainen & Luoto 2010, Kultti et al. 
2011). 

The sediment samples were prepared apply- 
ing standard methods for fossil midge analysis 
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described in Hofmann (1986), Walker (2001), 
and Brooks et al. (2007). The sediments were 
sieved through a 100 um sieve and the midge re- 
mains extracted from the residues using a 
Bogorov counting chamber (Gannon 1971) or a 
Petri dish under a stereomicroscope (25-40 times 
magnification). The fossil remains were perma- 
nently mounted in Euparal® or Canada balsam 
(Rhenohistol®) on microscope slides and identi- 
fied to the highest taxonomic resolution possible 
under a light microscope at 400-times magnifica- 
tion. A minimum of 100 midge individuals were 
identified from each sample, which exceeds the 
recommended minimum count size of 40-50 
head capsules for fossil chironomid analysis 
(Larocque 2001, Quinlan & Smol 2001). The 
identification was mainly based on the guides of 
Wiederholm (1983) and Brooks et al. (2007). 
Heiri et al. (2004) was used to identify the Tany- 
tarsini and Rieradevall and Brooks (2001) the 
Tanypodinae head capsules. The identification of 
Cricotopus pulchripes-type was based on the de- 
scription of Nyman et al. (2005). The phantom 
midges (Chaoboridae), biting midges (Ceratopo- 
gonidae), and black flies (Simuliidae) were iden- 
tified according to descriptions of Finnish speci- 
mens by Luoto and Nevalainen (2009), Luoto 
(2009b), and Luoto (2010). The black fungus 
gnats (Sciaridae) were identified according to 
Heiri and Lotter (2007). 

All statistical re-analyses of the previously 
available data were performed using relative taxa 
abundances and square-root transformation of as- 
semblage data to stabilize the variances of indi- 
vidual taxa expressed as percentages. Detrended 
correspondence analysis (DCA) was applied to 
choose between linear- or unimodal-based me- 
thods in further numerical analyses by estimating 
the lengths of compositional gradients (DCA axes 
1 and 2). The DCAs were run with detrending by 
segments and down weighting of rare species. 
DCA is an indirect ordination method that sum- 
marizes the variation in species assemblages 
along latent gradients (ter Braak & Smilauer 
2002). Canonical correspondence analysis 
(CCA) and redundancy analysis (RDA) were 
used to explore the relationships between midge 
assemblages and environmental variables. CCA 
is aunimodal ordination technique used with long 
gradient lengths and RDA is a linear ordination 


technique used with shorter gradient lengths. 
Both analyses can be used to identify environ- 
mental variables that are strongly related to the 
species assemblages (ter Braak & Smilauer 
2002). The CCAs and RDAs were run with only 
one environmental variable (without any trans- 
formations) at a time and by using statistically 
significant co-variables. The statistical signifi- 
cance of each variable was tested with a Monte 
Carlo permutation test (999 unrestricted permuta- 
tions), and variables were considered significant 
if the permutation test value (P) was < 0.01. In the 
intralake dataset, habitat (stream flow) and the 
presence of Nuphar lutea L. Sm., Equisetum 
fluviatile L. and Phragmites australis (Cav) Trin. 
ex Steud. stands were entered as nominal vari- 
ables, using binary codes (value 0 when absent 
and 1 when present). Habitat was entered as a 
“dummy” variable (0=lentic, 1=lotic), because 
discharge can vary greatly among or within sea- 
sons and following a precipitation event (Gan- 
doin et al. 2006). The stream flow analysis used 
all the 34 samples, but for other variables in the 
intralake dataset only the 29 lentic samples were 
used. The DCAs, CCAs and RDAs were per- 
formed using the program CANOCO, version 
4.52 (ter Braak & Smilauer 2002). 
Invertebrate-based cluster analysis (two-way 
indicator species analysis, TWINSPAN) was 
used as a divisive technique to group the lakes 
from the datasets and to identify indicator taxa for 
the environmental variables studied. TWIN- 
SPAN can be considered suitable for the present 
datasets, which consist of noisy data, many zero 
values, and high number of taxa (Gauch & 
Wittaker 1981). TWINSPAN was performed us- 
ing the program WinTWINS, version 2.3 (Centre 
for Ecology and Hydrology & University of 
South Bohemia, Huntingdon & Ceské Budejo- 
vicé) (Hill & Smilauer 2005). Spatial analyses 
were performed using Mantel non-parametric 
tests on the midge assemblages to determine if 
spatial autocorrelation existed in the intralake 
dataset, which may be considered as a strongly 
spatially structured environment. A Bray-Curtis 
distance was used for the biological distance ma- 
trix and geographical distance between the sites 
for the spatial distance matrix (1,000 randomi- 
sations), following the example of Quinlan et al. 
(2003). The geographical distances between the 
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Table 2. Summary of the statistically significant variables (P<0.01) affecting midge assemblages at different geo- 
graphical scales in Finland. The number of observation sites (N), results of the significance test (P), species-en- 
vironment correlations, and the proportion of variance explained. Species-environment correlation is the correla- 
tion between the site scores that are weighted averages of the species scores and the site scores and the site 
scores that are a linear combination of the environmental variables (Jongman et al. 1995). The regional, 
semiregional, and site-specific datasets were analyzed using the canonical correspondence analysis (CCA) and 
the local dataset using redundancy analysis (RDA). The values in brackets indicate results obtained by using 
statistically significant co-variables. In addition to the statistical tests, the gradient of the particular environmental 


variable in the dataset is shown. 








Scale Variable N P Correlation % variance Gradient 
explained 
Regional Mean air T 82 0.002 (0.001) 0.898 (0.905) 6.2 (6.4) 11.3-17.1 °C 
Water depth 82 0.002 (0.009) 0.750 (0.721) 2.4 (2.2) 0.5-9.0 m 
Semiregional DO 30 0.002 (0.005) 0.923 (0.881) 9.1 (6.2) 0.5-11.8 mg I” 
Conductivity 30 0.002 (0.040) 0.876 (0.868) 7.5 (5.2) 8-312 uS cm” 
Mean air T 30 0.002 (0.003) 0.817 (0.815) 6.7 (6.9) 15.9-17.1 °C 
Local Fish 9 0.001 (0.013) 0.969 (0.987) 37.0 (28) 0/1 FI* 
Water colour 9 0.004 (0.115) 0.900 (0.978) 29.2 (21) 25-125 Pt mg I’ 
Site-specific Stream flow 34 0.001 (0.001) 0.939 (0.913) 8.8 (7.8) 0/1 LI** 
Water depth 29 0.001 (0.001) 0.924 (0.920) 8.6 (7.2) 0.5-7.3 m 
Macrophytes 29 0.001 (0.082) 0.938 (0.900) 7.9 (4.6) 0-270 m (DLV***) 





*Fl = fish index, **LI = lotic index, ***DLV = sampling distance to littoral vegetation 


sampling points used coordinates, which were 
derived using a map tool because coordinates 
were not measured during the sediment sampling. 
Mantel test was performed using the program 
PAST (Hammer et al. 2001). 


4. Results and discussion 


The observed patterns of community structure 
can be affected by the scale of observation (Wiens 
1989, Levin 1992) or by the fact that some bio- 
logical phenomena are independent of scale 
(Aronson 1994). The ecological problems related 
to scale become more acute when studying tem- 
poral changes in community structure (Roem- 
mich & McGowan 1995) and their possible rela- 
tionship with anthropogenic effects (Smith & 
Buddemeier 1992). The temporal and spatial 
scales of observation are extremely important 
when interpreting changes in community struc- 
ture (Edmunds & Bruno 1996). 

The results from the dataset across Finland 
(regional) suggested that the mean air 7, , is an 
important factor affecting chironomid distribu- 
tion and abundance, together with sampling 
depth, whereas in the warm southern regions 


other factors also became significant (Table 2, 
Figs. 1, 2). The distribution of midge taxa in the 
regional dataset corresponded with general fau- 
nal observations made in Finland (Nyman et al. 
2005, Paasivirta 2009, Raunio et al. 2010). Of the 
individual taxa, Psectrocladius (Monopsectro- 
cladius) calcaratus-type, Micropsectra insigni- 
lobus-type, and Tanytarsus lugens-type preferred 
low temperature in the present dataset, while taxa 
such as Chironomus plumosus-type were more 
common in warm southern lakes (Fig. 2a). The 
temperature preferences were also evident in the 
TWINSPAN grouping of samples, which identi- 
fied indicator taxa representative for different 
parts of the climatic gradient in the regional 
dataset (Fig. 3a). Apparently, many of the taxa 
succeeding in particular limnological conditions, 
such as an elevated trophic status, may also fa- 
vour warm lakes (Fig. 3). It has been shown that 
in northern and alpine areas, nutrient conditions 
are related to summer temperatures, with cool 
lakes often having lower nutrient concentrations, 
while lowland areas at lower latitudes and alti- 
tudes are more often subject to anthropogenic nu- 
trient enrichment (Brodersen & Anderson 2002, 
Heiri & Lotter 2005). This is also probably one of 
the reasons for the variation found in the faunal 
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Fig. 2. Ordination plots of species (see Schnell et al. 

(1999) for standardized codes for taxon names) and 

the significant environmental variables using canoni- 

cal correspondence analysis. — a. Regional dataset. — 

b. Semiregional dataset. 


assemblages in the present study, as the warm 
southern parts of Finland are extensively utilized 
by humans, whereas the environment in the cold 
northern Finland is more close to its natural state. 
Unfortunately, the remote northern sites in the 
present study could not be examined for their wa- 
ter quality in as close detail as the southern lakes 
due to broken measuring equipment. 

In the dataset restricted to southern Finland 
(semiregional) with a temperature gradient of 
only 1.2 °C, parameters closely related to water 
quality, including dissolved oxygen and conduc- 
tivity, were the most important environmental 
factors associated with midge distribution and 
abundance (Table 2, Fig. 2b). Surprisingly, sum- 
mer air temperature was still a statistically signifi- 
cant correlate even with the very narrow tempera- 
ture gradient. Temperature, metabolism, and re- 
spiratory requirements are directly linked to each 
other, and respiration physiology is thus among 


the indirect causal responses and mechanisms 
making midges useful as indicators of climate 
(Brodersen et al. 2004, 2008). The oxygen condi- 
tions can change under stable temperature, but 
climate hardly ever changes without affecting the 
oxygen status (Brodersen & Anderson 2002). 
Under significant temperature changes, oxygen, 
productivity, and climate in combination are 
probably the forcing factors determining midge 
distribution. It has remained unclear whether chi- 
ronomids respond to temperature directly 
through physiological temperature tolerances or 
indirectly through changes in the lake habitat and 
other intralake processes (Anderson et al. 2008, 
Brodersen et al. 2008: Fig. 4). Nevertheless, pre- 
vious studies confirm the present results that lim- 
nological factors can have an overwhelming in- 
fluence in warmer climatic conditions or within 
short climatic gradients (Brodersen & Linde- 
gaard 1999, Brooks et al. 2001). However, of the 
limnological factors that were measured, pH was 
found not to have a significant influence on midge 
assemblages in the present regional and semi- 
regional datasets (Luoto 2009a, Luoto & Salonen 
2010). This is in contrast with the results of 
Nyman et al. (2005), who showed that pH was an 
important explanatory variable in their semi- 
regional dataset from north-western Finnish Lap- 
land. It should also be noted that we found no sig- 
nificant correlation between air temperature and 
measured water quality parameters in the semi- 
regional dataset. 

The present results indicate that climate may 
determine the midge assemblages at the regional 
scale, whereas water quality may become more 
important at the semiregional scale (Table 2, Fig. 
1). Interestingly, these parameters are closely re- 
lated as climate strongly regulates lake productiv- 
ity (Cohen 2003). Progressive eutrophication un- 
avoidably leads to changes in biota and biodiver- 
sity, blue-green algal blooms, hypolimnetic oxy- 
gen depletions, and fish kills. These changes ob- 
viously affect also the midge distribution through 
various feedback mechanisms at many levels. 
The potentially confounding phenomenon of 
midge trophic and climate optima being corre- 
lated may be alleviated when a dataset is properly 
interpreted in a scale-appropriate context. 

The TWINSPAN classification of the nine 
study lakes representing the local scale showed 
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Fig. 3. TWINSPAN grouping of lakes and the assigned indicator taxa (see Schnell et al. (1999) for standardized 
codes for taxa names) based on the midge assemblages in Finland. — a. Regional and semiregional datasets. — 


b. Local dataset. — c. Site-specific dataset. 


that the three lakes with no fish stocks were 
grouped together (Fig. 3b), suggesting that the 
faunal assemblages in fish-free lakes differ from 
those with fish, namely Perca fluviatilis, which is 
known to prey on chironomids (Diehl 1995). 
Most of the taxa were found in lakes with fish, but 
taxa such as Psectrocladius (P.) sordidellus-type, 
Psectrocladius (A.) flavus-type, Corynoneura 
lobata-type, and Chaoborus flavicans-type were 


completely restricted to fish-free sites. Chao- 
borus flavicans is known to co-occur with fish 
only in deep basins that provide refugia in the 
hypolimnion, in clay-turbid basins that prevent 
visual predation by fish, and in humic brown-wa- 
tered basins that have reduced fish predation po- 
tential due to their light-, oxygen-, and tempera- 
ture profiles (Liljendahl-Nurminen 2006). Inter- 
estingly, in addition to fish stock, water colour 
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was the only environmental factor that was found 
statistically significant in the present study (Table 
2) and the fish-free lakes in the dataset were those 
with dystrophic brown water (Nevalainen & 
Luoto 2010). However, when the CCA was run 
for water colour with the fish data as a co-vari- 
able, the influence of colour became statistically 
insignificant for midges. Because the number of 
samples in the local dataset was very low, the 
present results remain inconclusive and must be 
treated with caution. 

A Mantel test on the intralake dataset indi- 
cated that the midge assemblages exhibited statis- 
tically significant spatial autocorrelation (P = 
0.04), but this relationship was very weak (r = 
0.07). The results also showed that hydrological 
variables of stream flow and water depth were 
important factors controlling midge distribution 
when the study was scaled to a single lake and the 
aquatic vegetation patterns also contributed as a 
significant explanatory variable when the mar- 
ginal effects were examined (Table 2). The taxa 
richness in the surface sediments of Lake Pieni- 
Kauro and the Saavanjoki River (site-specific 
dataset) was very high, since over 100 midge taxa 
were encountered (Luoto 2010). This was most 
likely due to the great diversity of habitats avail- 
able for midges, including lentic and lotic envi- 
ronments with different exposures, substrata, and 
vegetation types. Furthermore, the midge assem- 
blages showed significant habitat-related hetero- 
geneity between the samples and the taxa were di- 
vided into clear indicator groups based on their 
habitat preferences in the intralake dataset (Fig. 
3c). The profundal zone in Lake Pieni-Kauro was 
dominated by the cold-indicating Heterotrisso- 
cladius maeaeri-type, while it was absent from 
other parts of the lake (Luoto 2010). Many other 
typical cold water taxa, such as Sergentia cora- 
cina-type, were also restricted to the profundal 
zone. The littoral areas of Lake Pieni-Kauro were 
characterized by taxa such as Cladotanytarsus 
mancus-type and Tanytarsus mendax-type, while 
the riverine samples were dominated by Thiene- 
mannimyia pseudocarnea-type and Simulium, 
fitting well with previous data on their ecology 
(Wiederholm 1983, Currie & Walker 1992, Gan- 
doin et al. 2006, Brooks et al. 2007). Psectro- 
cladius (P.) sordidellus-type and Cladopelma 
viridulum were strongly associated with Phrag- 


mites australis stands (Luoto 2010), which is in 
agreement with previous records of these taxa in 
Finland (Tolonen et al. 2001). 

The fact that some lakes in this study are in ur- 
ban areas and the semi-regional scale has more 
lakes affected by human activities than the re- 
gional scale does cause some bias due to uneven 
distribution of human impact. Furthermore, the 
environmental gradients are not the same be- 
tween scales and not all parameters could be 
tested between the datasets (e.g. climate vs. 
inlake processes), hence making the comparisons 
between them less successful. Although benthic 
communities can be sampled at almost any scale, 
not all environmental variables can be measured 
reasonably from all the scales. Nevertheless, the 
present implication that different causative pro- 
cesses are important at different scales is in agree- 
ment with the conclusions made by Walker et al. 
(1992). The ultimate environmental forcing fac- 
tors vary among different taxa, but it is also a mat- 
ter of spatial and environmental scaling; at the re- 
gional scale the ultimate forcing factor is proba- 
bly the prevailing summer temperature, at the 
semiregional scale water quality becomes more 
important, at the local scale fish predation can 
control the midges, and at the site-specific scale 
habitat characteristics determine the species com- 
position. However, many of these forcing factors 
are interconnected and the influence of a particu- 
lar parameter is often indirect and mediated 
through other environmental factors. It is also in- 
triguing that some of the parameters may be im- 
portant at several scales; summer air temperature 
at the regional and semiregional scales and water 
depth at the regional and site-specific scales. In 
the future, it would be important to aim to exam- 
ine midge distribution across spatial scales with 
all the parameters examined from all the scales. 

One ofthe most important implications in this 
study is that locally the predation pressure caused 
by fish may have overwhelming influence on the 
midge community composition. Therefore, an in- 
troduction or disappearance of fish, or just a 
change in the size of fish stock can potentially 
cause a very serious bias in any midge-based en- 
vironmental assessment. In addition to the influ- 
ence of food web changes, intralake processes, 
including changes in hydrology (water depth, 
stream flow) and in habitat availability (e.g. 
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macrophytes) may cause problems when examin- 
ing the determinants of a community and its dy- 
namics. At the different spatial scales, commu- 
nity structure is driven by several interconnected 
factors such as local history, different rates of 
change of community structure and the patchy ef- 
fects of stochastic events (e.g. disturbances and 
recruitment). However, in accordance with pre- 
vious studies (e.g. Sorebón & Nakamura 2009), 
abiotic conditions are fundamental at large scales 
and bionomic effects at smaller scales that makes 
a theoretic background to link the different de- 
scriptions of species distribution across scales 
and the factors influencing them (Hortal et al. 
2010). The scenopoetic factors, which become 
neglible at site scales, are responsible of number 
of processes influencing species distributions, in- 
cluding physiological constraints, responses to 
climate and habitat gradients, active habitat selec- 
tion, and range shifts in response to change in cli- 
mate or habitat (Hortal et al. 2010). Furthermore, 
previous studies have confirmed that climate can 
regulate the ranges of species at large scales, 
while variables related to habitat operate at local 
and site-specific scales (Pearson & Dawson 
2003, Thuiller et al. 2004). To understand species 
distributions and the effects of forcing factors, it 
is necessary to consider the effect of scale (Hui et 
al. 2010). It is also important to understand that 
no single scale can be considered “correct” for 
system analysis (Levin 1992) and extreme cau- 
tion should be used when scaling the results to 
scales that were not sampled. 
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